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Abstract 

We present results of electronic structure calculations for well-relaxed Au/benzene- 
1,4-dithiol/Au molecular contacts, based on density functional theory and the gen- 
eralized gradient approximation. Electronic states in the vicinity of the Fermi energy 
are mainly of Au 5d and S 3p symmetry, whereas contributions of C 2p states are 
very small. Hybridization between C 2p orbitals within the benzene substructure is 
strongly suppressed due to S-C bonding. In agreement with experimental findings, 
this corresponds to a significantly reduced conductance of the molecular contact. 

Key words: molecular contact, electronic structure, electrical transport 
PACS: 73.20.-r, 73.20.At, 73.40.-c, 73.40.Cg, 73.63.Rt, 85.65. +h 



1 Introduction 



Molecular devices nowadays are becoming more and more important because 
of a rapidly increasing number of potential applications, spanning a wide range 
from new computer architectures to multifunctional chemical sensors and med- 
ical diagnostic tools. Exploring the use of specific molecules as active compo- 
nents in electronic circuits has been a field of special interest in recent years. 
Electrical properties of molecular devices are usually addressed by measuring 
the transport through a molecule sandwiched between metallic electrodes. The 
resulting transport characteristics, of course, reflects both the metal-molecule 
contact and the molecule itself. Because transport in general is restricted to a 
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small number of valence orbitals, the features of a molecular device strongly 
depend on the local electronic structure at the different molecular sites. De- 
tails of the chemical bonding, such as hybridization or charge transfer, hence 
play a key role for understanding the device, and structural relaxation due to 
the metal-molecule contact has to be taken into consideration. 

Experimentally, Au/benzene-l,4-dithiol/Au transport interfaces are prepared 
via the mechanically controllable break junction technique [TJ. In a first step 
a thin gold wire, covered with a self-assembled benzene-l,4-dithiol monolayer, 
is streched until breakage. The opposing gold contacts then are slowly moved 
together until the onset of conductance is reached. From the theoretical point 
of view, the prototypical Au/benzene-l,4-dithiol/Au interface likewise has at- 
tracted much attention. Band structure methods have succeeded in reproduc- 
ing the shape of the experimental current- voltage characteristic of the contact 
when an external bias is applied. However, the absolute value of the conduc- 
tance is overestimated by one or two orders of magnitude [2.3.4J. The current 
magnitude appears to be an extremely sensitive function of both the contact 
geometry and chemistry. But even if atomic details of the contact are taken 
into account, transport calculations fail to reproduce the experimental fact of 
a conductance gap at zero bias. 

The benzene- 1,4-dithiol molecule has ir bonding and n* antibonding orbitals 
formed by the carbon and sulfur p orbitals perpendicular to the ring plane, as 
well as a bonding and a* antibonding orbitals due to in-plane orbital overlap. 
While both the a and 7r orbitals are fully occupied, the antibonding orbitals 
stay empty. Because they mediate direct orbital overlap, a states are subject 
to stronger bonding- antibonding splitting than n states, and therefore appear 
at lower energies. For a small bias, the conductance of the molecular contact 
shows Ohmic behaviour because of a smooth density of states near the Fermi 
energy. When the bias is increased, two conduction peaks are observed, which 
can be traced back to resonant tunneling through either ir* or (at higher bias) 
ir orbitals [2]. 

The erroneous prediction of metallic transport characteristics for experimen- 
tally insulating molecules has been associated with the continuous exchange- 
correlation approximation applied in standard density-functional schemes [5]. 
For the Au/benzene-l,4-dithiol/Au contact, it has been demonstrated that an 
atomic self-interaction correction is useful to compensate this shortcoming in 
an approximate way, opening a conduction gap in the current-voltage charac- 
teristic. However, in this letter we show that the modification of the chemical 
bonding in the benzene ring, resulting from the specific contact geometry, has 
serious effects on the local electronic structure and thus likewise explains the 
insulating state. We thus conclude that taking into account the very details of 
the crystal and electronic structure is mandatory in order to obtain adequate 
input for transport calculations. 
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2 Calculational Details 



The electronic structure calculations presented in the following are based on 
density functional theory and the generalized gradient approximation. To be 
specific, we use the WIEN2k program package, a state-of-the-art full-potential 
code applying a mixed lapw and apw+lo basis [6J. In particular, the WIEN2k 
program allows for structural optimization of the interface geometry. Thus it 
is suitable to study effects of covalent bonding and hybridization in detail. A 
well-relaxed interface geometry is of great importance for reliable band struc- 
ture results, as even small alterations of the orbital overlap, due to structural 
relaxation effects, may cause significant modifications of the electronic states 



In our calculations, the charge density is represented via 88491 plane waves, 
the k-mesh for the Brillouin zone integrations comprises 72 k-points, and the 
Perdew-Burke-Ernzernhof parametrization of the exchange-correlation poten- 
tial is applied. While Au 5p, S 3s, and C 2s orbitals are treated as semi-core 
states, the valence states comprise Au 5d, 6s, dp, S 3p, 4s, C 2p, 3s as well 
as H Is orbitals. Corresponding atomic sphere radii amount to 2.6 ae for the 
Au, 1.8 as for the S, 1.3 for the C, and 0.7 ag for the H spheres. 

\ y 




Fig. 1. Schematic structure of the Au/benzene-l,4-dithiol/Au molecular contact. 
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Table 1 

Bond lengths in the relaxed Au/benzene-l,4-dithiol/Au molecular contact. 

Our calculation relies on the contact geometry shown in figure [TJ The benzene- 
1,4-dithiol molecule is sandwiched between two gold contacts in fee [001] ori- 
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Fig. 2. Au/benzene-l,4-dithiol/Au molecular contact: Partial Au and S densities of 
states. The DOS is normalized with respect to the number of contributing sites. 

entation. While the carbon site C2 forms the center of the molecule, the site 
CI contacts the sulfur atom on top of a gold pyramide. Each sulfur atom thus 
has four nearest neighbour gold sites. Moreover, the fourth gold layer off the 
contact forms the basal plane of our unit cell, which allows us to use periodic 
boundary conditions. In total, each electrode consequently comprises 23 gold 
atoms in four layers. A convenient starting point for the contact geometry is 
given by the bulk gold bond lengths and angles, the bond lengths of the bare 
molecule, and a typical Au-S distance. Structural relaxation is carried out by 
means of minimization of interatomic forces as well as the total energy. The 
resulting bond lengths for the optimized contact geometry, as summarized in 
table [U largely resemble the numbers reported by Basch and Ratner [12]. Our 
subsequent discussion of the Au/benzene-l,4-dithiol/Au molecular contact is 
based on electronic structure data for the fully relaxed geometry. 



3 Results and Discussion 

Results of our LAPW calculation for the Au/benzene-l,4-dithiol/Au interface 
are shown in figures [2] and EJ We present partial densities of states (DOS) for 
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Fig. 3. Au/benzene-l,4-dithiol/Au molecular contact: Partial C densities of states 
for carbon sites 1 and 2 (see figure [Q). The DOS is normalized with respect to the 
number of contributing sites. 

the contact Au sites (forming the top of the gold electrodes), the S sites, and 
the two carbon sites CI and C2; compare the structure sketch in figure HJ For 
a quantitative analysis, the density of states is normalized with respect to the 
number of contributing sites in each case. In the energy range shown, almost 
all electronic states trace back to Au 5d, S 3p, and C 2p orbitals, respectively. 
The partial DOS of the contact Au sites is dominated by a broad structure of 
Au 5d states, which is centered at roughly — 3.5eV and extends far above the 
Fermi energy. The gross features of the DOS shape resemble the findings for 
fee gold, reflecting metallic conductivity for these sites. Moreover, hardly any 
deviation from fee gold is found for the electronic structure of Au sites in the 
third gold layer off the contact, which justifies the setup of our electrodes. 

The S 3p DOS is rather unstructured in a wide energy window reaching from 
about — 5eV to 4eV. Especially, at the Fermi level a distinct DOS is present, 
with a magnitude similar to the Au 5d DOS. This observation traces back to 
hybridization between gold and sulfur orbitals, as expected for Au-S bonding. 

In the next step we turn to the CI site, which is affected by S-Cl and C1-C2 
bonding. As compared to the previous findings for gold and sulfur, we have a 
significantly reduced (but still finite) DOS in the vicinity of the Fermi energy. 
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Moreover, strong hybridization between CI 2p and S 3p orbitals is evident, as 
both densities of states have sharp peaks near — 9eV and — 8eV, and common 
structures around both — 6eV and — 4eV. Hence, remarkable S-Cl bonding is 
a typical feature of the benzene- 1,4-dithiol molecule. 

Sharp DOS peaks at —9 eV and —8 eV and a distinct structure reaching from 
about — 7eV to — 5eV likewise are present for the C2 site. These states belong 
to a bonding orbitals of p symmetry, which is also true for the corresponding 
S and CI states. Of course, a bonding orbitals mediate the main part of the 
orbital overlap within the benzene ring and with adjacent sulfur atoms. They 
therefore are subject to strong hybridization. States at energies above — 4eV 
are expected to belong to tt bonding and n* antibonding orbitals. In fact, the 
C2 2p DOS has characteristic peaks near — 3eV and 2.5 eV, which, however, 
are not found in the CI 2p DOS. Hybridization between n bonding orbitals at 
sites CI and C2 hence is almost negligible. As a consequence, 7r-type coupling 
within the benzene ring is seriously reduced in case of the benzene- 1,4-dithiol 
molecule, which we attribute to strong S-Cl bonding. Due to the latter, CI n 
and Ti* states are subject to energetical shifts, yielding a rather unstructured 
DOS at energies above — 4eV. Since tt bonding states are responsible for the 
electrical conductance of the benzene ring, reduced C1-C2 coupling leaves the 
whole ring insulating. In fact, the C2 2p DOS reveals hardly any states in the 
vicinity of the Fermi energy. 

In conlusion, we have studied the electronic structure of the Au/benzene-1,4- 
dithiol/Au transport interface by means of band structure calculations within 
density functional theory. Taking into account the specific interface geometry, 
we have addressed the chemical bonding in the molecular contact. While cou- 
pling via a bonding S 3p and C 2p orbitals is strong, hybridization between n 
bonding orbitals at neighbouring carbon sites is almost completely suppressed 
due to S-C bonding. As a consequence, the local DOS at carbon sites without 
sulfur bonding almost vanishes at the Fermi energy, which explains the exper- 
imental observation of an insulating state. A transport calculation neglecting 
the details of the chemical bonding in the benzene-l,4-dithiol molecule is ex- 
pected to overestimate the conductance (at least at small bias), although the 
shape of the current- voltage characteristic may be predicted correctly [2|3|4] . 
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